Sterile neutrino decay as a common origin for 
LSND/MiniBooNe and T2K excess events 
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We point out that the excess of electron-hke neutrino events recently observed by the T2K collab- 
oration may have a common origin with the similar excess events previously reported by the LSND 
and MiniBooNE experiments and interpreted as a signal from the radiative decays of a sterile neu- 
trino Vh with the mass around 50 MeV produced in u^ neutral current (NC) interactions. In this 
work we assumed that the vh can also be produced in i/tNC reactions. 

PACS numbers: 14.80.-j, 12.60.-i, 13.20.Cz, 13.35.Hb 



Over the past 10 years there is a puzzle of the 3.8 
a event excess observed by the LSND collaboration [l|. 
This excess originally interpreted as a signal from F^ -^ 
77g oscillations was not confirmed by further measure- 
ments from the similar KARMEN experiment J2[. The 
MiniBooNE experiment, designed to examine the LSND 
effect, also found no evidence for z^^ — >■ v^ oscillations. 
However, an anomalous excess of low energy electron- 
like (e-like) events in quasi-elastic neutrino events has 
been observed [3| . New MiniBooNE results from a search 
for 77^ — ^ Fg oscillations also show an excess of events, 
which has a small probability to be identified as the 
background-only events [4|. The data are found to be 
consistent with V^ -^ Fg oscillations in the 0.1 eV^ range 
and with the evidence for antineutrino oscillations from 
the LSND. 

In the recent work [5| (see also PtS]) it has been shown 
that puzzhng LSND, KARMEN and MiniBooNE results 
could all be explained in a consistent way by assuming 
the existence of a heavy sterile neutrinos {I'h)- The Vh 
is created in u^^ neutral-current interactions and decay 
subsequently into a photon and a lighter neutrino v in 
the LSND and MiniBooNE detectors, but it cannot be 
produced in the KARMEN experiment due to the high 
energy threshold. The v^ could be Dirac or Majorana 
type, and it could be produced, e.g. through the v^ — Vh 
mixing. The Vh could decay dominantly into jv pair if, 
for example, there is a large enough transition magnetic 
moment between the Vh and v nrass states. Such kind 
of VhS nray be present in many interesting extensions of 
the standard model, see e.g. 9j. Assuming the Vh is pro- 
duced through mixing with f^, the combined analysis of 
the LSND and MiniBooNe excess events suggests that the 
Uh mass, mixing strength, and lifetime are, respectively, 
in the range 



40 < m,, < 80 MeV, lO'^ < \Uf,h\'^ < 10^ 

10"" <Th< 10^^ 



(1) 



A detailed discussion of consistency of these values with 
the constraints from previous searches for heavy neutri- 
nos [lO| as well as of the interpretation of the Vh decay 
in terms of transition magnetic moment is presented in 
[3, Q . Briefly, the mixing of ([T]) is not constrained by the 



limits from the most sensitive experiments searched for 
extra peaks in two-body tt, K decays [10], because the i^h 
mass range of H]) is outside of the kinematical limits for 
7r^2 decays, and not accessible to K^2 experiments due 
to experimental resolutions. The parameter space of ([T]) 
cannot be ruled out by the results of high energy neu- 
trino experiments, such as NuTeV or CHARM, as they 
searched for i/h 's of higher masses {rrih > 200 MeV) de- 
caying into muonic final states {fiTTiy, iifj-v, fiev,..) [10l |. 
which are not allowed in our case. The best limits on 
\U^h\^ derived for the mass range H]) from the search for 
Vh — > e^e~v decays in the PS191 experiment [Ilj, as well 
as the LEP bounds |12| , are found to be compatible with 
P^ assuming the dominance of the Vh decay. New limits 
on mixing \U^h\^ obtained by using the recent results on 
precision measurements of the muon Michel parameters 
by the TWIST experiment [IJl are also found to be con- 
sistent with ([TJ. Finally, the most stringent bounds on 
|C/^/iP coming from the primordial nucleosynthesis and 
SN1987A considerations, as well as the limits from the 
atmospheric neutrino experiments, are also evaded due 
to the short Uh lifetime. 

Very recently, the T2K collaboration, which study u^^ 
neutrino neutrino interactions in a long baseline experi- 
ment at J-PARK, has reported on observation of an ex- 
cess of electron- like events in charge-current quasi-elastic 
(CCQE) neutrino events over the expected standard neu- 
trino interaction events |1J| . A confirmation of the T2K 
excess and clarification of its origin have great impor- 
tance for neutrino physics. Although the most popular 
mechanism for this excess is v^ — ^ v^ neutrino oscilla- 
tions with nonzero value of the neutrino mixing angle 
013, one can still reasonable ask if neutrino oscillations 
are the only explanation for the T2K result, see e.g [l5|. 
In this work we study a possible manifestation of the 
presence of VhS in the J- PARC neutrino beam and show 
that the excess of e-like events observed by T2K could be 
interpreted as a signal from the production and radiative 
decay of a Uh previously suggested for the explanation of 
the origin of similar excess events observed by the LSND 
and MiniBooNe experiments. 

In the T2K experiment, specifically designed to search 
for Vn — >■ i/g oscillations and for measurements at the first 



t'/i — Vt oscillation maximum (corresponding to the atmo- 
spheric neutrino parameters sin^2023 = 1 and Am^^ — 
2.4 X 10"'^ eV^), the neutrino flux at the far detector lo- 
cation is dominated by j^t's. Therefore, we will assume in 
the following that the i^h can be also produced in z^,-NC 
interactions, e.g. through the mixing with the Vt neu- 
trino [16] . Taking into account that the corresponding 
mixing strength \Urh\'^ is purely constrained by exist- 
ing experimental data, makes this assumption more in- 
teresting. For example, the bounds from the NOMAD 
experiment |C/rftp < (4 — 1) x 10~^ for the v^ masses 
from 40 to 80 MeV were obtained from the search for 
the Vh — i> VtG'^g~ decay under assumption that this de- 
cay is the principal decay mode of the Vh- These limits 
can be significantly relaxed assuming the dominance of 
the Vh — >■ ^v decay. Direct searches for radiative decays 
of an excited neutrino v* — >■ v^ produced hi Z —^ u*!/ 
decays have also been performed at LEP [l^l. For the 
best limit Br(Z -^ v*v)By{v* -^ 1/7) < 2.7 x lO"*^ from 



ALEPH [12], taking into account g^t^- 
we find [3| 



IC/. 



\Urh\^ X 



^ [MeV] 



[■s] 



< 4.8 X 10^ 



(2) 



For the mass and mixing range of ([T]) it results in r^,^ > 
10"" s. 

In addition, there is also a hint from the measurements 
of the D+ -^ T+Vr decay rate at CLEO [H [Ilj. It has 
been observed that for some decay modes of outgoing 
T-leptons, the obtained decay rate of -D+ — s> r'^Vr is sig- 
nificantly higher than the predicted one, so that com- 
bined branching ratio deviates from the theoretical pre- 
diction at the level of ~ 10%. One may speculate that 
this inconsistency is due to existence of a sterile neutrino 
mixed into the tau neutrino. This would result in an 
additional contribution to the £>+ -^ t^v^ decay rate 
from the decay D+ -^ T^Vh- It has been found that 
for the mass range below 190 MeV the mixing strength 
required to explain D^ — >■ t'^Vt discrepancy should be 
\Urh\'^ = 0.16 ± 0.09 0, which for the mass and lifetime 
range of ([T]) is consistent with ^ within the large error. 

The T2K experiment is described in details in [20]. It 
uses an almost pure off-axis v^ beam originated from the 
7r+ and K decays in flight, which are generated by 30 
GeV protons from the J-PARC Main Ring accelerator. 
The detector consists of a near detector complex, used 
to measure precisely the v^ flux and to predict the stan- 
dard neutrino interaction rate in the far detector, which 
is the Super-Kamiokande (SK) water Cherenkov detector 
located at the distance of 295 km from the proton tar- 
get. The SK detector is a cylindrical tank, about 40 m in 
diameter and 40 m height, filled with ~50 kt of purified 
water 21!]. The detector has a fiducial volume (FV) of 
22.5 kton within its cylindrical inner detector (ID), and 
a 2 m-wide outer detector (OD) served as a veto against 
cosmic rays and neutrino interactions in the surround- 
ing rock. The Cherenkov light rings generated by muon. 



electron and converted photon tracks are used for the 
reconstruction of the events. The T2K search for e-like 
events from t'p — > z^e neutrino oscillations uses the data 
sample collected during the years 2010-2011 14]. The 
strategy of the analysis is to identify the i^eCCQE can- 
didate events by reconstructing in the FV isolated single 
e-like rings that are accompanied by no other activity 
in the outer detector. The measured rate of the e-like 
events is then compared to the one expected from known 
reactions. An excess of AiV = 4.5 electron-like events (6 
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FIG. 1: Schematic illustration of the production and subse- 
quent radiative decay of a heavy neutrino in the SK detector 
(top view). The VhS are produced in fpf^-jNC interactions 
of the J-PARC Uf_t neutrino beam either in the FV, the ID 
region outside FV, the OD or, as shown, in the surrounding 
rock. In the later case the v^ would penetrate the rock shield 
and would be observed in the SK neutrino detector through 
their Uh — >■ "yv decays followed by the decay photon conversion 
into an e^e~ pair in the ID FV. 

events observed and 1.5±0.3 expected) has been observed 
in the data accumulated with 1.43 x 10^" protons on tar- 
get (pot). For the following discussion several distinctive 
features of the excess events are of importance [ij: a) 
the excess is observed for single e-like tracks, originating 
cither from an electron, or from a photon converted into a 
e^e^ pair with a typical opening angle ~ me/i?e+e- < 1 
degree (for E^+e- > 100 MeV), which is too small to 
be resolved into two separate Cherenkov rings in the SK 
detector(here, me,ii^e+e- Eire the electron mass and the 
e+e^ pair energy); b) the reconstructed neutrino energy 
is in the range 200 < E^^ < 1000 MeV. The vari- 
able E^^ is calculated under the assumption that the 
observed electron track originates from a i^eCCQE in- 
teraction; c) the visible energy E.^is is required to be 
Evis ^ 100 MeV; d) the angular distribution of the ex- 
cess events with respect to the incident neutrino direc- 
tion is wide and consistent with the shape expected from 
i/gCCQE interactions; e) there is no additional significant 
activity in the OD detector. 

To satisfy the criteria a)-e), we propose that the ex- 
cess events are originated from the production and sub- 
sequent radiative decay of a heavy neutrino Vh in the 
FV of the SK detector. The heavy neutrinos are as- 



sumed to be produced in the neutral-current quasi-elasic 
(NCQE) interactions '^^{r) + iV — >• i^^ + iV of muon or 
tau neutrinos either in the SK FV, the ID region outside 
FV, OD region, or in the surrounding rock. In the later 
case, if v^ is a relatively long-lived particle, the flux of 
i/h's would penetrate the rock shielding without signifi- 
cant attenuation and would be observed in SK through 
their Vh -^ 1^ decays with the subsequent conversion 
7 -^ e+e~ of the decay photons in the SK water target, 
as schematically illustrated in Fig. [TJ Similar to v^ — ^ v^, 
neutrino oscillations, the occurrence of vy^ — > 71/ decays 
would appear as an excess of single e-like events from de- 
cay photon conversion in the SK detector, above those 
expected from standard neutrino interactions. To make 
a quantitative estimate, we performed simplified simula- 
tions of the Vh production and decay processes in the SK 
discussed above. 

The flux of the produced VhS can be calculated by 
using the following equation for the Vh production cross 
section, 

(3) 
where a{v^(^T-) + A^ — > ^'^(r) + N) is the cross section for 
i^^(t-)NCQE interactions and / is the phase space factor 
calculated for this two-body reaction, which takes into 
account dependence on the Vh mass. The energy spectra 
of the produced i^/i's, whose momenta pointing to the SK 
fiducial volume, as well as the angular distribution of the 
I'/i's, were calculated for different Vh masses by taking 
into account the v^, v^ energy distributions at the far de- 
tector [l4J. In these simulations we used a parametrized 
i^^ energy spectrum obtained at far detector from the 
reconstructed i^^tCCQE events [14]. The Vr energy dis- 
tribution at the FD position was calculated from the pri- 
mary v^ spectrum within the standard two-neutrino os- 
cillations scheme for the atmospheric parameters quoted 
above. The total number of NC events in the FV of the 
SK detector was used for normalization. 

Once the Vh flux was known, the next step was to cal- 
culate the e+e" spectrum based on the v^ — >■ ^v decay 
rate. For a given flux ^^^ , the expected number of signal 
events from v^ — > ^v decays occurring within the fiducial 
length L of the SK detector with the FV entrance point 
located at a distance L' from the v^ production vertex is 
given by 



E 



j,,d^.. 



dEi^i^ ) v^^{v^) 



exp( ^ 



(4) 



[1 _ exp(-:^;^^)] ^e^s,+,-dE,^dV 

Pvh'^h i-tot 

where pu^ is the Vh momentum and th is its lifetime at the 
rest frame, Fg+g- , Ttot are the partial and total mass de- 
pendent t-zi-decay widths, respectively, e is the e+e~ pair 
reconstruction efficiency and the integral is taken over 
the FV, ID region outside FV, OD and the surrounding 
rock volume. It is assumed that the total rate Ttot of the 
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FIG. 2: The number of i^h decays from the Vh's produced in 
the rock as a function of rock thickness calculated for i/,t and 
i^T energy spectra at the far detector, assuming \U^h\ ~ 1 
and If/rhP = 1, for different th values shown in the plot, 
a = 0, and 1.43 x 10'^'' pot. The fraction of i^h ^^ 7^^ decays 
from z^^NC reactions is < 15%, while the rest is due to z/tNC 
interactions. 



i^h decays is dominated by the radiative decay Vh -^ 7^^, 
i.e. the branching fraction BR{vh -^ 7^^) = ('^fa-»7t^) ^ ^ 
Q . The acceptance A of the SK detector was calculated 
by tracing the produced VhS to the detector FV taking 
momentum and angular distributions into account. The 
energy of the photon from the Vh decay depends on the 
initial neutrino energy and on the center-of-mass angle O 
between the photon momentum and the Vh momentum 
direction. Therefore, the photon laboratory energy spec- 
trum depends on the cm. angular distribution, which 
is generally given by dN/dcosQ ^ 1 + a ■ cos8, where 
asymmetry coefficient a is in the range —1 < a < 1 for 
Dirac, and a = for Majorana neutrinos |22l |. The recon- 
struction efficiency of the photon converted in the fiducial 
volume of the SK detector was taken to be ~ 70% from 
the T2K simulations of the i/g QECC events [13] • An 
example of the calculated number of the Vh decays in the 
SK FV is shown in Fig. [2] as a function of the thickness 
of the rock surrounging the detector for the i^h mass of 
50 MeV and several Th values. It is seen, that if the Vh is 
a relatively short-lived particle, i.e — — < 1 the number 
of the signal events is quickly saturated and does not de- 
pend on the rock thickness. Neutrino interactions, with 
little hadronic activity in the final state, occuring in the 
OD or SK support structure, as well as in the part of the 
ID outside FV, can also yield an isolated e-track from the 
i^h decay in the FV. The attenuation of the Vh - flux due 
to Vh interactions in the rock with the average density 
3.2 g/cm^ was found to be neghgible. 

In Fig. [3] an example of distributions of the kinematic 
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FIG. 3: Distributions of the excess events reconstructed as 
j^eCCQE events in the SK detector as a function of variable 
E^^ for Evis > 100 MeV from the experimental data sample 
(dots), and from a combination of the Vh — >■ 7^^ decay of z^h's 
produced in i^rNC plus expected neutrino background (bot- 
tom shaded histogram, from Ref. [1J|) calculated for a = — 1 
(solid line), a = (dashed line), and a = +1 (dashed-dotted 
line) cases shown in the plot, the fh mass of 50 MeV, and the 
fh lifetime Tj,^ = 10~^ s. Error bars include only statistical 
errors. The distributions are normalized to six events, which 
corresponds |C/rhP — 0.025. A distribution of neutrino back- 
ground plus u^ —>■ Ue neutrino oscillations at sin^20i3 ~ 0.1 
(dotted histogram, from Ref.[ij]) is also shown for compari- 
son. 



variable E^^ for the excess fh decays events in the SK 
detector reconstructed as jvCCQE events plus neutrino 
background predicted in [ij] are shown for Eyis > 100 
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MeV, m„ 



50 MeV, 



10 s, and different values by n^^ ~ nfj,\U^h\ 



of a. These distributions are calculated for the domi- 
nant production of v^s by z^'^-'s assuming that the e~^e~ 
pair from the converted photon is mis-reconstructed as 
a single track from the z^eCCQE reaction. The distri- 
butions are then normalized to six events, which corre- 



sponds lUrh 



0.025, to compare them with the T2K 



data. Simulations are in reasonable agreement with the 
experimental distributions. For instance, for the distri- 
butions shown in Fig. [3] the x^ test of their consistency 
with T2K data yields p - values of 0.79, 0.87, and 0.92 
for a = —1, 0, and -|-1, respectively. In these calcula- 
tions only statistical errors for both experimental, as re- 
ported in 14], and simulated spectra are included. The 
test is based on the method of comparison of experi- 
mental and simulated histograms recently proposed in 
[23| , which can be apply for analyzing data samples con- 
taining small numbers of events per bin. The simulated 
excess events, shown in Fig. |31 are mainly distributed 
in the region 200 < E^^ < 1200 MeV. The simula- 
tions showed that the shape of the E^^ distributions 
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FIG. 4: Number of expected e-like events from Uh decays 
plus neutrino background in the SK FV for 1.43 x 10^° pot as 
a function of the mixing strength It/rhP and the Vh lifetime 
Th from the region of ([T]) calculated for ruh = 50 MeV and 
a = 0. 



is sensitive to the choice of the i^h mass, t^,^, and, in 
particular, the a-parameter values: the smaller the a, 
the softer the spectrum. The distribution of cosine of 
the opening angle between the e-like ring and neutrino 
beam direction is found to be consistent with j/g CCQE 
events. Taking into account the estimated number of 71 
events expected to be observed in the FV from the ordi- 
nary NC interactions [1J|, the total number of Vh events 
inside the fiducial volume of the SK detector is given 
It^-hP, where coefficients n^^rir 
- 70 - 160, de- 
pending on Th- As |C/^iip < It/T/il"^ by a factor of a few, 
mainly i^t-NC interactions in the FV and in the rock con- 
tribute to the total number of excess events. The number 
of excess events from interactions in the FV is, roughly, 
oc It/r/iP, while the number of events from the rock is 
oc \Urh\'^Th- In Fig. [4] an example of distribution for 
the expected number of signal plus neutrino background 
events is shown in more details in the {tu] \Urh\'^) parame- 
ter space. For the larger mixing the number of Vh events 
increases, while for smaller lifetime values it decreases 
due to the more rapid decays of i^/i's, and mostly NC in- 
teractions in the FV contribute in this case. The fraction 
of excess events expected to be seen in the OD depends on 
the Vh lifetime and typically is < 20% (< 1 event) which 
is consistent with observations [ij]. It is also interesting 
to compare the spatial distributions of the excess events, 
which, in the case of Vh decays, is a combination of the 
uniform distribution from the NC interactions in the FV 
and a distribution from NC interactions outside the FV. 



vary in the range n^ ~ 12 — 30 and 
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FIG. 5: The probability to observe more than 5 events in 
the T2K experiment as a function of ^^l.h\^ and IC/rhl^ for 
Th = 10-^ s. 



The later is the z/^ hfetime dependent: for shorter t^ the 
excess events are expected to be distributed presumably 
near the edge of the FV. One interesting feature of the 
spatial distribution of events in SK is that they are all in 
the first half of the detector. To reproduce this feature 
one could try to use the Vh lifetime shorter than t^ < 0.5 
ns, so that i/^'s produced in the rock decay mainly in 
the first part of the SK. For shorter lifetimes, as one can 
see from Fig. [21 the number of events from the rock de- 
creases, and mixing values |C/t/i|^ ^ 0.03 are required to 
keep the signal in the range of 4-5 events. Note, that 
the increase of If/r/iP results also in increase of the rel- 



ative fraction of events produced in NC interactions in 
the FV (typically, about 40% for the long r/,.'s), which 
are distributed uniformly over the FV. 

The probability P{ny^ > 5) to observe more than 5 
events in the T2K experiment calculated as a function 
of mixings |C/pftp and \Urh\'^ is shown in Fig. [5] It was 
obtained by using the approach of Ref.[24] and taking 
into account the uncertainties in the background esti- 
mate [lj|. For example, for P{ni,^ > 5) > 0.25, using 
^ to constrain the Th, the most suitable values of the 
parameters are 

40 < m^ < 80 McV, 10"^ < \Uf,h\^ < 10"^ 



10"^ < \Urh\^ <3x 10" 



10" 



<T„<10-9s. (5) 



In summary, in this work we study a possible manifes- 
tation of the presence of heavy neutrinos in the J- PARC 
neutrino beam and show that, assuming the Vh mixing 
into the i/t, the T2K excess events could originate from 
the same mechanism as those observed by the LSND and 
MiniBooNE experiments, namely from the production 
and radiative decay of a sterile neutrino with properties 
of ^ . This interpretation is found to be compatible with 
all the constraints a)-e). The distribution of the excess 
events in kinematic variable E^^ is found to be consis- 
tent with of the shapes of distributions obtained within 
this interpretation. Our analysis may be improved by 
more accurate and detailed simulations of the T2K ex- 
periment, which are beyond the scope of this work. A 
definite conclusion on the presence of Vh events can be 
drawn when the T2K statistics is substantially increased. 
Finally, note that several ideas on searching for Vh in /z 
decays [8j, with existing neutrino data [25|, in radiative 
K decays [2^, [23| , and with neutrino telescopes [2g have 
been recently proposed. The author thanks D.S. Gor- 
bunov, N.V. Krasnikov and M.E. Shaposhnikov for use- 
ful discussions and/or comments, and S.I. Bityukov, A.E. 
Korneev and D. Sillou for help in calculations. 
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